Mesoscopic Behavior Near a Two-Dimensional Metal-Insulator Transition 
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We study conductance fluctuations in a two-dimensional electron g function of chemical 

potential (or gate voltage) from the strongly insulating to the metallic regime. Power spectra of the 
fluctuations decay with two distinct exponents {1/vi and 1/vh). For conductivity a ~ 0.1 e'^/h, we 
find a third exponent (l/«i) in the shortest samples, and non-monotonic dependence of Vi and vi 
on a. We study the dependence of Vi, vi, Vh, and the variances of corresponding fluctuations on a, 
sample size, and temperature. The anomalies near a ~ 0.1 e^ /h indicate that the dielectric response 
and screening length are critically behaved, i. e. that Coulomb correlations dominate the physics. 
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The metal-insulator transition (MIT) is one of the fun- 
damental problems in condensed matter science. Re- 
cent theoretical work strongly suggests the crucial 
role played by electron-electron interactions in the tran- 
sition regime. However, since the development of the 
scaling theory of localization it has been asserted 
that all states are localized in 2D, in agreement with 
early experiments |^ on relatively low- mobility samples. 
A recent experiment on a two-dimensional electron 
gas (2DEG) in Si metal-oxide-semiconductor field-effect 
transistors (MOSFETs) provides evidence for the exis- 
tence of a true MIT. The samples used in that exper- 
iment [Q had a much higher mobility and, as a result, 
the Coulomb interactions played a greater role relative 
to disorder. Thus it has been speculated Q that this 
MIT is driven by interaction effects. Using mesoscopic 
measurements, we provide direct evidence for the crucial 
role of Coulomb interactions at the MIT in a 2DEG. 

We investigate the statistics of conductance fluctu- 
ations in a 2DEG as it undergoes a transition from 
strongly insulating to metallic behavior. In the insulat- 
ing regime, electrons move in a strong, random potential 
by tunneling through localized states. In our relatively 
small samples and at low temperatures, the total number 
of states that contribute to conduction can be small, e. g. 
of the order of 50-100. By sweeping the gate voltage Vg, 
the chemical potential is shifted relative to the energy 
of localized states. As a result, the number and nature of 
localized states that dominate the transport change, and 
the conductance G changes up to several orders of magni- 
tude. In addition, it is well established that, deep in the 
insulating regime, the density of states D{E) increases 
exponentially with increasing Vg |^ . 

Our measurements were carried out on n-channel 
MOSFETs fabricated on the (100) surface of Si doped 
at « 3 X 10^'' acceptors/cm'^ with 500 A gate oxide 
thickness and (unintentional) oxide charge < 10^° cm~^. 
The peak mobilities of our samples were of the order of 
2 m^/Vs - comparable to those exhibiting a true MIT [Q . 



In contrast to those samples, ours are much smaller: rect- 
angular with sourcc-to-drain lengths L — 1 — 8 /im, 
and widths W = 11.5 — 162 /im. They were short 
enough to exhibit conductance fluctuations, and wide 
enough to ensure good statistics. Measurements were 
performed in a dilution refrigerator with a lock-in at 
^ 10 Hz and an excitation voltage of 0.2 /iV rms. Con- 
ductance was measured as a function of Vg at tempera- 
tures 0.01 < T < 0.8 K. 

The behavior of the typical conductance, i. e. of 
(InG) (averaging over intervals of V^), in these samples 
has been studied in detail ||^. At very low Vg and for 
T < 0.1 K, the conduction is due to resonant tunnel- 
ing through channels containing a few localized states. 
At higher T, the conduction proceeds via variable-range 
hopping along isolated chains of several hops, consistent 
with the model in Ref. . The conductance of each chan- 
nel depends not only on the energy and position of each 
localized state in the channel but also on the Coulomb 
interactions between electrons. At higher Vg, the current 
paths become more complicated as more and more states 
contribute to conduction. 

The conductance fluctuates on two, and sometimes 
three different scales oi Vg These scales are ap- 

parent in the raw data (see Fig. 1(a) inset), but a de- 
tailed statistical analysis is necessary for inferring the 
underlying physics. We have analyzed the power spectra 
of fluctuations 51nG = InG(Vg) - (lnG(K,)) for differ- 
ent Vg ranges or conductivities cr = (L/W) exp(lnG), 
sample sizes, and temperatures. The power spectrum 
is a Fourier transform of the autocorrelation function 
C{AVg) ^ {S\nG(yg)d\nG(yg + AVg)). A typical power 
spectrum S{l/AVg) is presented in Fig. |l](a). There 
are two distinct regions of exponential decay, S{f) = 
5(0) exp(— 27rw/), each corresponding to a Lorentzian 



G{AVg) with characteristic widths 



The voltage 



correlation scales vi and Vh represent the average peak 
spacings of the corresponding fluctuations. They are ob- 
tained by fitting both the low- frequency (/ < /o) and 
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FIG. 1. Power spectrum S{f) of fluctuations S\nG. (a) 
L = 2.0 fj.m, W = 162 ^m; T = 0.023 K. The vertical dashed 
line separates the low-frequency and high-frequency regions. 
The fluctuations 51n G vs. Vg, corresponding to this S{f), are 
shown in the inset, a = 1.37 x 10"^e^//i. (b) T = 0.045 K. 
The vertical dashed line separates the low-frequency region 
from the rest of the spectrum at approximately the same fre- 
quency /o as in (a). A new energy scale emerges at inter- 
mediate frequencies (located between the two vertical lines in 
the plot) in the transition regime between strongly localized 
and metallic behavior. The correlation voltage at the highest 
frequencies remains unchanged. 



high-frequency (/o < /) parts of S{f) separately by the 
above exponential form, /g, which separates the low- 
frequency part from the rest of the spectrum, does not 
depend on any of the parameters within the scatter of 
our data; vq — l//o — (0.45 ± 0.15) mV, averaged over 
all of our measurements. Fig. [|(b) illustrates the evolu- 
tion of the spectrum as a increases. While the correlation 
voltage at the highest frequencies remains essentially un- 
changed (~ 50 /iV) over a wide range of a, we observe the 
emergence of a new distinct correlation voltage Vi (i. e. a 
new Lorentzian in C{AVg)) at intermediate frequencies. 
From Fig. ^ which plots the location of the boundary 
fc between the intermediate-frequency and the highest- 
frequency regions (shown by the dash-dotted vertical line 
in Fig. |l](b)), we see that this new feature of S{f) is most 
pronounced for 10~^ < cr (e^//i) < 10""'^. It vanishes 
when the system becomes either more insulating or more 
metallic. 

^From the work on quantum dots and wires and 
tunneling through localized states ||], it is known that 
Coulomb interactions give rise to peaks in conductance 
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FIG. 2. Characteristic voltage scale Vc = 1/ fc, correspond- 
ing to the boundary /c between the intermediate-frequency 
and highest-frequency regions of the spectrum, as a function 
of a. The dashed line is a guide to the eye. 
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FIG. 3. High-frequency (/o < /) correlation voltage vs. 
a for samples of different size. Vi (/o < / < fc) is shown 
with solid symbols, and Vh {fc < f) with open symbols. One 
sample (A) had W — 162 /im, and the rest had W = 11.5 /im. 
Dashed lines guide the eye. 



such that the peak spacing is determined by the strength 
of the interaction. Therefore, at least some of the mea- 
sured correlation voltages should be proportional to the 
typical Coulomb energy in the system. As the MIT is ap- 
proached from the insulating side, the Coulomb energy 
U in our samples increases as a result of a decrease in 
the mean separation of the electrons. As the MIT is ap- 
proached from the metallic side, U also increases because 
the screening becomes less efficient P, p^ , pT| . Therefore, 
we expect U to have a maximum in the transition region. 

Obviously, it is important to analyze the evolution of 
the correlation voltages with cr in more detail. We start 
with the high-frequency (/o < /) part. Fig. ^ presents 
Vh and Vi for several samples of different size. Fits to the 
decay rate are shown for fo < f < fc (solid symbols) and 
for fc < f (open symbols). For fc < f, this correlation 
voltage Vh remains of the order of ^ 50 /iV, independent 
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FIG. 4. Variance C(0) of high-frequency (/o < /) fluctua- 
tions vs. (J at several temperatures. Solid symbols are from 
the fo < f < fc part of the spectrum, and open symbols from 
the fc<f region. W — 11.5 ^m. 
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FIG. 5. Correlation voltage i;; of low-frequency (/ < /o) 
fluctuations vs. a for samples of different size. One sample 
(A) had W = 162 y,m, and the rest had W = 11.5 /xm. The 
lines are guides to the eye. 



of sample size, T, and even g. For the shortest samples, 
however, the third correlation scale Vi (solid symbols in 
Figs. H and ^) appears near a ^ 0.1 e'^/h. Vi is observ- 
able only in the smallest samples, and it also disappears 
at higher T, so it seems that this feature of S{f) van- 
ishes as a result of averaging. Vi goes through a dramatic 
(increasing by 400%) maximum near a = 0.1 e^//i. 

In order to characterize the fluctuations completely, we 
have also analyzed their variances C(AVg = 0). Fig. ^ 
shows how C(0) of the fluctuations for /o < / depends on 
a in our smallest sample. For both types of fluctuations 
{fo < f < fc and /c < /) in the high-frequency range, we 
observe a power-law decrease of C (0) with an increasing 
a. For the highest-frequency fluctuations, we flnd no T- 
dependence of the variance in any of our samples but 
we have evidence that C(0) is reduced by increasing the 
sample area. Although the variance of the intermediate- 
frequency fluctuations seems to decay with increasing T, 
we do not have sufficient data to determine the functional 
form. For this part of the spectrum, C(0) seems to be 
suppressed rapidly by an increase of sample area. 

In the low- frequency range (/ < /o), the correlation 
scale vi has a maximum as a function of a (Fig. ||), in a 
qualitative agreement with the expected behavior of U. 
The maximum occurs again for 10~^ < a (e^/ft.) < 10^^. 
It is suppressed by increasing the sample size. The re- 
sults obtained on two L = 2 fim samples with differ- 
ent W point to the absence of area averaging in vi. A 
rise in temperature shifts the peak to higher values of 
a. For example, in our smallest sample {L = 1.0 /xm and 
W = 11.5 ^m) the peak moves from a « 2.5x 10~^ (e^//i) 
at T = 0.010 K to cr « 2.0 X 10~^ (e^/h) at T = 0.555 K. 
It should also be noted that vi is about 1-2 orders of 
magnitude larger than Vh- C'(O) for low- frequency fluc- 
tuations does not depend on any of the parameters that 
we have varied except that there seems to be a small 
decrease of C(0) as a increases beyond 10^^ {e^ /h). 



The rate of the chemical potential shift with Vg is ||] 
dn/dVg = {dfi/dn,){dn,/dVg) = {l/D{E)){C/e), where 
Us is the 2D electron density, and C is the gate capaci- 
tance of our devices (6.92 x 10~^ F/m^). Since it is not 
possible to relate D{E) to a without an adequate theo- 
retical model of transport, we are not able to obtain, for 
example, an accurate dependence of various correlation 
energies on a. It is possible, however, to make some esti- 
mates, assuming that D{E) varies from its constant 2D 
value of 1.6 X 10^^ cm^^eV^^ in the metallic regime down 
to about 10-20% of that value in the insulating limit. 

With this in mind, we find that the voltage correlation 
scale of ~ 50 /xV, corresponding to the highest-frequency 
regime, varies from about 1.35 fieV at low a, to 0.15 fxeV 
at the highest a. In spite of a relatively weak dependence 
of correlation energy on a, these fluctuations disappear 
eventually at high Vg because of a decrease in their am- 
plitude. It is interesting that these fluctuations are well 
resolved all the way up to 0.8 K even though their cor- 
relation energy is much smaller than ksT which, while 
not intuitively plausible at flrst glance, is quite reason- 
able upon detailed consideration of this issue for a re- 
lated problem |lj]. Given that this is the smallest en- 
ergy scale, hardly dependent on any of the parameters, 
and that these fluctuations are so ubiquitous, we are led 
to speculate that their correlation energy is related to 
the energy level spacing in our system. If true, it would 
correspond to a sample area of about 2x2 ^m. This is 
somewhat smaller than the area of our devices but it is 
consistent with the fact that only a fraction of the sample 
dominates conduction, as discussed in the introduction. 
We recognize that there are other energy scales of impor- 
tance in this problem, but without a detailed model for 
the transport, we must forgo further speculation. 

In the same fashion, we find that the correlation en- 
ergy of intermediate-frequency fluctuations, which ap- 
pear only in the crossover regime, is of the order of 0.5- 
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1 /xeV, also smaller than fc^T. For low-frequency fluctu- 
ations, the correlation energy varies from about 25 /zeV 
in the strongly insulating regime, to 45 /xeV in the tran- 
sition region (assuming D{E) ~ 10^''cm~^eV~"'^), down 
to 5 /xeV at the highest a. Therefore, it appears that it 
has a maximum in the crossover regime. It is comparable 
to or larger than ksT. In addition to having obtained 
qualitative agreement with the expected behavior of a 
Coulomb interaction energy, we can also make a reliable 
quantitative comparison in the most insulating regime, 
where U in our devices is reduced only by screening by 
the metallic gate Q. Then U ^ 25 ^eV corresponds to 
a reasonable distance of 0.3 /im, comparable to the mean 
separation ('^ 0.8 /zm) of localized states lying within an 
energy of a typical level spacing of /i. Using the same 
expression for U in the case of intermediate-frequency 
fluctuations, J7 ~ 1 /ueV gives a similar distance of 1 /im. 
Even though, strictly speaking, this estimate is not quite 
accurate in the transition region, it could indicate that 
these fluctuations reflect some more subtle properties of 
the Coulomb interactions very close to the MIT, - per- 
haps related to the beginning of overlap between chains 
of states or more metallic islands - as a result of the 
competition of several important length scales. 

In addition to three distinct correlation voltages, 
we have identified two characteristic voltage scales, 
vq and Vc- Vc is a measure of the width of the 
intermediate-frequency regime, observable most clearly 
in the crossover region. In energy units, Vc varies from 
about 15 fj,eV at low a, to 0.6 fieV in the transition re- 
gion, up to about 1.5 /ieV at high a. This indicates the 
existence of a minimum of this characteristic energy as 
a is varied. Our current data also show a tendency for 
the minimal value of Vc in Fig. ^ to go to zero as T is 
lowered but more detailed measurements of the temper- 
ature dependence would be useful. On general grounds, 
one expects that the characteristic length scales, such as 
the localization length, will diverge at the MIT, and that 
the associated characteristic energy will go to zero. Our 
results for Vc thus strongly support the existence of a 
MIT in this 2D system. 

In summary, we have analyzed the power spectra of 
fluctuations (5 In G as a function of Vg. None of the non- 
interacting theories, which might be applicable to our 
experiment fl^ within a limited range of a, provide an 
adequate description of our flndings, either because of the 
failure to predict the correct shape of the low-frequency 

or because of the failure to 



part of the spectrum |15 
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describe the observed temperature dependence of the cor- 
relation voltages 0jl^. But most importantly, general 
considerations based on non-interacting models sug- 
gest monotonic behavior of all characteristic energy scales 
as a function of n^. In this paper, we present a striking 
flnding of a strong nonmonotonic behavior of several en- 
ergy scales, which may be understood by considering the 
effects of Coulomb interactions. Our experiment makes 



it absolutely clear that a theory that treats both disorder 
and electron-electron interactions on an equal footing is 
required in order to describe the MIT in this 2D system. 
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